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Abstract 

We consider extensions of excluded volume interactions for com- 
plex corpora that generalize simple rod-like particles. The Onsager 
equation can be defined for quite general configuration spaces, and 
the dimension reduction of the phase space in the limit of highly in- 
tense interaction can be shown. The formalism describes both freely 
articulated and interacting N-rods and the example of interacting 2- 
rods is given in detail. 



1 Introduction 

When rod-like particles are suspended in fluids, the particles can be described 
by relatively simple configuration spaces M, for instance, the unit sphere 
M = S 1 in R 2 or M = §> 2 in M. 3 . The probability of finding a rod whose 
director m belongs to the region EcMis J s f(m)dm. The measure dm on 
the unit sphere is the uniform measure, the natural volume element (length 
or area) induced from the ambient Euclidean space. The probability measure 
f(m)dm characterizes the particle distribution. In equilibrium, this measure 
is obtained by minimizing the free energy (|15j) 

S\f\ = J [log f + ^U^f dm 
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where b is a nonnegative parameter, representing a combination of inverse 
temperature and intensity of interaction, and 



is a potential computed from an interaction kernel u(m,p). The physical 
modeling is reflected in the choice of kernel. Only two effects are taken into 
account in the free energy: an entropic effect, representing thermal fluctua- 
tion of the particles, and a quadratic mean field effect representing excluded 
volume interactions between particles. The description of the problem is 
completed by specifying the configuration space M, the uniform measure dm 
and the interaction kernel u. Once these are determined we may ask ques- 
tions, such as: do minima of the free energy exist, how many such minima 
exist, what happens to them as the parameter b varies from b = to infinity. 

In the case of more complicated corpora, the configuration space M can 
be rather complicated. We use the word "corpus" to refer to a body that 
has finitely many degrees of freedom, such as an assembly of articulated 
rods, or rods connected to balls of different sizes and with different motion 
constraints between the parts. The collection of all corpora of a particular 
problem is the configuration space. It is useful to phrase the equilibrium and 
kinetic problems broadly, in quite general configuration spaces. We consider a 
metric space M, with distance d, and assume we are given a Borel probability 
measure p on M. This probability measure corresponds to the normalized 
uniform distribution of corpora m G M, and M is the configuration space. 
For absolutely continuous probability Borel measures v « p, the potential 
function 



associated to v is defined by a kernel u. The function u(x,y) is given by 
the model, and it is real valued, symmetric in x, y, Lipschitz continuous and 
bounded. Let f(x) be the Radon-Nikodym density / = ^ so that 





v = fdp,, 



f > 0, p— a.e., f M fdp = 1, and let associate to it a free energy 




(1) 
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Because f M Ufdfi = f M J M u(x,y)f(x)f(y)dfj,(x)dLi(y), there is no loss of 
generality by assuming that u(x, y) is symmetric, as only the symmetric part 
\{u(x, y) + u(y, x)) contributes to the free energy Critical points of the free 
energy obey 

58 

which leads to Onsager's equation: 

/ = Z- l e~ w (2) 

Indeed, this follows taking the Gateaux derivative (i.e. the first variation) of 
£ in a direction h, i.e. 

^-£[f + th] lt=0 = 0, 

with arbitrary h, subject only to J M hdp, = in order to respect f M (f + 
th)dfi = 1. A simple calculation, using the fact that u(x,y) is symmetric in 
x, y leads to J M (log / + bUf)hdii = for all h. This implies that log f + bUf 
is constant /z-almost everywhere. Because u is Lipschitz continuous, it is easy 
to see that the solutions of ([2]) are Lipschitz continuous, positive and never 
vanish. 

The kernel u gives the rule to construct the potential U associated to the 
particle distribution v = fd/x; in examples this is computed from physics 
and it instructs the particles to conform to each other. As b — > oo there 
are phase transitions. In the abstract case where we have a distance d we 
may consider general u(x,y) = $>(d(x,y),x,y) and the simplest examples 

are u(x,y) = d(x } y) p with p = 1,2 or u(x,y) = —e ~~P . The interaction 
potential u(x, y) should have a minimum at x — y. 



2 Rods 

The simplest example of single rods has M = S 1 . The Onsager kernel is 
u(9i, 62) = I sin(#i — 2 )\, and the Maier-Saupe kernel is 



2/ 



u(9i, 6 2 ) = — cos 

Kernels are defined up to addition of constants, so the Maier-Saupe kernel 
can be taken as u(6i, 6 2 ) = - c ° e t 2 ( 6l -^i)) _ 
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Throughout this paper we will use the notation JCf = —U[f] = —U to 
emphasize the linear dependence of U on /. (The minus sign is a convention 
we maintain to be consistent with previously published literature.) 

In the Maier-Saupe potential case, JCf projects / on the eigenspace of 
the Laplace-Beltrami operator corresponding to eigenvalue 4. If the function 
/ is written as 



1 1 

f(e) = - + -J2y^os(2 3 e) (3) 

/TT TT ' * 

then Maier-Saupe potential is 



2n n 

3=1 



JCf =i|frcos(20). (4) 
The solutions of the Onsager equation in that case are of the form 

g (r) = (Z(r))- 1 e ra »W (5) 

with 

i-2-k 

Z{r) = / e rcos{2e) d9 (6) 
Jo 

and in Fourier representation 

1 1 °° 

<?(r)(0) = - + cos(2j0) (7) 

3=1 

Onsager's equation is equivalent to the implicit transcendental equation 

2r 

9i(r)) = J- (8) 

If b < 4 this equation has one solution, namely r = 0. If b > 4 there is 
exactly one more solution for r > at r = r(b). There is exactly one solution 
for r < 0, at r = — r(b) corresponding to a rotation of | of the solution with 
r > 0. It is known also that go(r) = 1 for all r, 

gj(r) > 0, fl- i+ i(r) < g^r), j > 1, r > 0. (9) 

If r = r(6) is determined by (jSJ), then 

fli(r(6)) = 1 - ^ 
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there are recursion relations to compute all gj(r), and 




holds for the inverse function, b(r). This is an increasing unbounded function 
for b > 4. The stability of the solution g(r) also follows from the inequality 
9i ~ 92 > 0. Both the inequality and the stability are not obvious. The limit 
of zero temperature is a delta function: 

lim g( r ( b )){9)d0 = 5 . (11) 

b — >oo 

The choice of 9 = is dictated by the symmetry 9 \— > —9 that we imposed. 
Otherwise, we obtain any delta function on the circle. The case of S 2 with 
Maier-Saupe potential is similar. The potential projects on the eigenspace of 
eigenvalue 6 of the Laplace-Beltrami operator on the sphere. Two implicit 
transcendental equations determine the solution. There is a finite tempera- 
ture phase transition and the zero temperature limits are delta functions on 
points or geodesies. 



3 Kinetics 

The kinetic description of rod-like particles ([9]) can be naturally generalized 
when the configuration space M is a Riemannian manifold. When M is 
Riemannian manifold, the kinetic equations are 

d t f = A 9 f-b&iv g {fV g {K,f)) (12) 

with A g , divg, V ' g Laplace-Beltrami, divergence and gradient in M, JCf given 
by 

JCf = - u(m,p)f(p)dfx(p) (13) 

J M 

and the uniform measure dfi is the Riemannian volume element. Note that 
the equation can be written as 

$/ = div p (/V,(log /-&£/)) 

The solutions are smooth, positive and normalized so they have unit integral. 
The free energy 
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is a Lyapunov functional: 

J t S = -J /|V g (log/-WC/)| 2 d,i(p) 

M 

If M is connected, the only possible steady solutions are solutions of On- 
sager's equation 

/ = Z~ x e hK K (15) 

The dynamical system is dissipative: the solutions are bounded after an 
initial transient time. The bounds, in very strong norms, are independent of 
the initial data. The global attractor is compact, finite dimensional and is 
formed with solutions of Onsager's equations and their unstable manifolds. 
In the case of M = S 1 with Maier-Saupe potential, the kinetic equations are 
given by the sequence of ODEs 

j t Vi = -^fyj + bjyi (y^i - y j+1 ) . (16) 

The equations on § 2 are more complicated. Some of the results concerning 
equilibria and kinetics for the Maier-Saupe potential can be found in pQ -[I], 
[B], [7], [TO] -[II] and [HI]. Recently it was shown that the system has inertial 
manifolds in both S 1 and § 2 ([UJ-[T8]). 

4 Freely articulated corpora 

We consider corpora made of N articulated rods that are allowed to rotate 
freely. The configuration space is M = §" x • ■ • x §" with n = 1 or 2. The 
potential is a sum of binary interactions 

N 

v(pi,qi,Pi,q2,'-PN,QN) = ^2 u j(Pj^j)- ( 17 ) 

3=1 

The uniform measure \l on M is the product measure. In this situation, the 
corresponding solution of the Onsager equation is a product measure 

u(d Pl . . . dp N ) = I^ =1 Zj 1 e~ bU ^dp j (18) 

where each Z~ 1 e~ bU ^ p ^ is a solution of Onsager's equation in the j compo- 
nent, and everything can be reduced to the study of individual rods. This 
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situation holds in general, when the corpora are made up of iV components, 
themselves corpora of a simpler kind, each with configuration space Mj and 
uniform measure \iy The composite corpora belong to the configuration 
space M = Mi x ■ ■ • x Mat, and if the kernel has the form ffTTl) then the 
corresponding operator 

fcfipi, ■ --Pn) = - I v (pi,gi, . . . ,PN,qN)f(qi, ■ ■ - ,gjv)^(gi) • --dii(q N ) 

J M 

is a sum of corpus-by-corpus operators, corresponding to like-parts only: 

N N 

icf( Pl , ... PN ) = Y, (**/) (Pi) = - E u Af\iPj) 

3=1 3=1 

with 

fcjfiPj) = ~ J u j(Pj' Qj)f(Qu ■ ■ ■ gjv)nd//j(ft). 

In this case, if we consider any solution of Onsager's equation, / = Z~ 1 e bK f , 
then the partition function 

Z= [e bRI Hdfii 
Jm 

is factored 

Z = Uf =1 Zj, with Zj = [ e b ^ f dfi( Pj ) 

and consequently the solution of Onsager's equation is a product of corpus- 
by-corpus solutions, 

f(Pi, ■■■Pn) = /i(pi)/2(pa) • • • Jn(pn) 

with fj(pj) = Z~ 1 e~ bUj( - p ^ a solution of the Onsager equation on Mj with 
uniform measure fij and interaction kernel Uj. The measure corresponding 
to the collective is the product measure of the single corpus solutions: the 
components interact only with like-components. In this case the study of 
the collective behavior reduces to the study of single corpus statistics, and 
in particular, the zero temperature limit can be any combination of zero 
temperature limits, component-wise. 
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5 A simple example: two interacting rods 



Clearly, if the components of complex corpora do not interact only with 
like-components, then a product measure is not the solution of the Onsager 
equation. Let us start by investigating a simple example. Assume that two 
rods with configuration space M = S 1 are articulated, so together the corpora 
belong to M = S 1 x S 1 . We will assume that the interaction between pairs 
of two-rods is determined entirely by the area of the triangle formed by each 
two-rod, and each two-rod corpus rejects two-rods that have very different 
area than themselves. The simplest interaction that achieves this is 

w(Pi,9i,P2,92) = \\e(pi) A e(p 2 ) - e(q 1 ) A e(q 2 )\\ 2 

with e(p) = (cos p, ship) if p e [0, 2n]. Each S 1 is viewed as a subset of 
IR 2 ; the exterior product /\ 2 R 2 is isomorphic to R. Every element in it 
is a multiple of e\ A e 2 , with e\ = (1,0) and e 2 = (0,1). For instance, 
e(pi) A e{p 2 ) = sin(p 2 — Pi)( e i A e 2 ). The norm satisfies 

||e(pi) A e(p 2 ) - e(qi) A e(g 2 )|| 2 = (sin(pi - p 2 ) - sin(gi - q 2 )f 

and consequently, the integral operator is 

Kfipufr) = - I \\e(pi) Ae(p 2 ) - e(qi) A e(q 2 )\\ 2 f(q u q 2 )dq!dq 2 . 

J M 

We take the uniform measure on S 1 x S 1 to be the natural measure dq\dq 2 . 
We note that the potential has the form 

U(pi,p 2 ) = sin 2 (j9i -p 2 ) - 2^sin(pi - p 2 ) + 7 

with z, 7 constants determined by /. Onsager's equation / = Z~ 1 e blc f re- 
duces therefore to 

/ z = [sin0](z,7) 
\ 7=[sin 2 0](^, 7 ) 

where we use the notation 

M(z, 7 )= I** <K9)g{9)M 
Jo 

q(Q\ — ^-lg-bsin 2 (6»)+26z sin 6-67 
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2tt 

—b sin 2 (6)+2bz sin S-b-y^g 

The solution is f(pi,P2) = g(pi ~ P2)- Note that g does not depend on 7. 
Let then 

s(0, z) = sin 6> — 

and let 

\ s ](b z) = Jo f^lfjf ™ 

The Onsager equation is equivalent to 

[s](6,2)=0. (19) 

This determines 2, which in turn determines g, f. Note that z = always a 
solution that yields 

/ ( Pl ,p 2 ) = Z- 1 e- fesin2 ^- p2 ). 

As b — » 00 this solution tends to — P2) mod7r), a degenerated two-rod. 
Consider now 

/•27T 

A(z,r) =65 / e- b{sine ~ z)2 d6 
Jo 

with r = b~ l . Note that 
so f|T9|) is equivalent to 

<9 2 A = 0. (20) 
Note also that A solves a linear heat equation with temperature as time: 

d T \ = \dl\ 

4 2 

The function A is even in z, so it is enough to study it on [0, 1]. The "initial 
value" obeys, for any z G [0, 1), 

lim X(z, r) = 2\/n- 
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This is an increasing function for z G (0, 1), the derivative is positive. Because 
the derivative obeys the same linear heat equation, and the heat equation 
preserves positivity, we might have been tempted to think that z = is the 
only solution of d z X = 0. But the behavior is more subtle because, clearly, by 
direct inspection, ||(l,r) < 0! The reason for this is the singular behavior 
at z — 1. As r — ► or b —>■ oo, we have a change of behavior of the function 
in a small region of the order 6~ 1//2 near z — 1. In this region there is a 
transition to much higher values of A, and the derivative d z X changes sign. 
This implies that there is a phase transition at positive r, i.e., there exists 
< z{b) satisfying 

d z X(z b ,r) = 
and lim T ^ z b = 1- Consequently 



is a delta function concentrated on right two-rods that make a right angle. It 
is instructive to note that the reason for this transition is that the vanishing 
of the gradient of the phase sin# — z can occur, when z — 1 at a minimum 
of the potential. Further details will presented elsewhere ([!]). 

6 A few general observations 

More complicated configuration spaces arise when the corpora are n-gons in 
space or in the plane. The natural conformation distance between such cor- 
pora is Hausdorff distance, modulo rotations. This area is rather open for 
investigations, and it is useful to start by stating the main general expecta- 
tions. It is expected that, if the configuration spaces are compact (a realistic 
assumption if the corpora have finitely many degrees of freedom and finite 
extensivity) , then generically, the zero temperature limit will be a singular 
measure, i.e., it will be concentrated on a set of zero /i measure, where /i is 
the uniform measure. In many examples the set of zero measure is just one 
point in M. 

Let us consider a general compact metric space M with distance d and 
"uniform" distribution /i, a Borel probability on M. Let us assume that 
there exist 0<&<l,c>0 




fi(B(x, r)) > ce 



(21) 
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for all x G M, and all r sufficiently small. Here B(x,r) is the ball centered 
at x G M and of radius r in the metric d. This assumption says that 
all balls are charged at least a small positive amount. Note that if M is 
any compact Riemannian manifold of dimension n and \x is the Riemannian 
volume element, then the condition (J2TJ) is automatically satisfied, because r n 
is much larger, for small r, than the right hand side of fl2~Tj) . The interaction 
kernel we consider is a function of the distance, u(d(x,y)), so 

U(x)= [ u(d(x,y))f(y)dn(y). (22) 

J M 

We assume that u is non-negative, bounded and Lipschitz continuous, i.e. 
there exist positive constants C and L so that 

< u{d) < C (23) 

and 

\u(d 1 )-u(d 2 )\<L\d 1 -d 2 \ (24) 
hold for all d,di,d 2 > 0. We assume also that 

u(0) = 0. (25) 

As we mentioned before, the interaction kernels are defined up to additive 
constants: if we add c to the interaction kernel, then the potential is changed 
by the same amount c and the free energy is changed by adding ^; its critical 
points, and in particular its minima are unchanged. Let us take a ball of B 
of radius r, set \ = [i(B)~ 1 1b the normalized indicator function of B, and 
compute 

£[ X ] = logGu(S)- 1 ) + \v(B)- 2 [ [ u(d(x,y))dfi(x)dfx{y) 



B J B 



We obtain, using (T2T| 1231123)1 that 



£\X}< ( I) +log( - c )+bLr 



which implies, by choosing r = b , that 



inf£[/] < b k + d 
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holds for large b with some constant C\. This means that uniform measure, 
whose energy is linear in b, 

= o / / u(d(x,y))dfi(x)dfi(y) 

1 J M J M 

does not achieve the minimum of energy for large b. On the other hand, 
for may examples, the uniform measure is a solution of Onsager's equation. 
Indeed, if T is a fi measure preserving isometry of M, then the potential 
associated to the uniform measure n, U[l] = —/CI, is invariant under right 
composition with T, i.e., 

U(Tx) = U(x) 

holds for all x G M and all measure preserving isometries. If measure pre- 
serving isometries act transitively, i.e., for any x,y G M there exists T a 
measure preserving isometry that maps x to y, y — Tx then, U = —JC(1) is 
a constant. This implies that the uniform measure is a solution of Onsager's 
equation. This is the case for many homogeneous spaces. The combination of 
these two very simple observations leads to the conclusion that a phase tran- 
sition occurs, whenever the homogeneous measure is a solution of Onsager's 
equation, the condition (1211) holds and the interaction kernel is a normalized 
Lipschitz function of distance. That simply means that, while for b = ob- 
viously the only solution of Onsager's equation is the homogeneous measure, 
and while this continues to be a solution for b > 0, at large enough b there 
exist other solutions as well, in quite great generality. 

Now we describe a general tendency of solutions to concentrate. Let 
dv = fdfx be any probability measure absolutely continuous with respect to 
/i and let U be the potential associated to it via (|22|) . Then, in view of the 
property (12"51) we have 

< U{x) < C (26) 
and from (j24|) and the triangle inequality we deduce that 

\U(x)-U(y)\<Ld(x,y) (27) 

so that the potentials associated to any probability fdfi are non-negative, 
uniformly bounded and Lipschitz continuous. 

Theorem 1 Let M be a compact metric space with distance d. Let fi be a 
Borel probability measure on M that satisfies (Tflj). Let u satisfy l[2^24^ . 
Then: 



12 



(A) For any b > there exists a solution g that minimizes the enery: 

£[g] = min £[f] 
f>o, f M fdn=l 

The function g solves the Onsager equation 

g(x) = (Z(6))- 1 e- 6l/ ^ 

with 

Z(b) = [ e- bU ^d/j(x) 

JM 

and 

U(x) = / u(d(x,y))g(y)dn(y) 

JM 

The function g is normalized f gdfi = 1, strictly positive and Lipschitz con- 
tinuous. 

(B) Let b n —>■ oo and let dv n = g n dfi be a sequence of solutions of Onsager 
equations corresponding to b n . By passing to a subsequence we may assume 
that the sequence converges weakly to a probability measure v = lim n u n . 
There exists a non-negative Lipschitz continuous function U^x) on M such 
that v is concentrated on the set 

E = {x e M I f/oofx) = miriyzMUooiy)} 

Thus, for any continuous function <fi supported in the open set M \ S, 

lim / (J)(x)g n (x)dfi = 



M 



Proof. For the proof of (A) we fix b > and note that S[f] is bounded below 
uniformly for all / > 0, j M fdfi = 1. We take then a minimizing sequence 

fr 

a = inf £[f] = lim S[fA. 

Without loss of generality, by passing to a subsequence and relabelling, we 
may assume that the measures fjdfi converge weakly to a measure dv. Using 
dSHlEZl) and the Arzela-Ascoli theorem, we may pass to a subsequence, which 
we relabel again fj , so that Uj converge uniformly to a non-negative Lipschitz 
continuous function U. Then it follows that 

U(x)= / u(d(x,y))dv(y) 

JM 
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holds and 



lim / Ujfjd/j,— / Udv. 
Jm Jm 

Because £[fj] is a convergent sequence of numbers, it follows that 



lim / fj log fjdfi 

exists. In particular, the above integrals are bounded uniformly. It then 
follows that dv is absolutely continuous i.e., dv = gdfi, with g > 0, g e 
L^idfx). Indeed, the sequence fjd/j, is uniformly absolutely continuous. This 
follows from the convexity of the function y log y and the Jensen inequality 

(n(A))- 1 / /log/d/x > ralogm 

J A 

where m = fi(A)^ 1 f A jd\i. Then we have 

mlogm < C/fj,(A) 

with a fixed constant C > 0, uniformly for all / = fj and any A. Let us 
choose R so that RlogR = C/fi(A). If we denote by / = j A fdfi, then either 
m < R, or, if not, then mlogi? < C//i(A). In either case the inequalities 
imply 

I<C/\og(R) 
and as fi(A) — > 0, R — > oo. This inequality signifies 



with lim^^o^l^) = and S(x) independent of n. This implies that v is ab- 
solutely continuous, v = gdfi with < g G The weak convergence 
tested on the function 1 implies that J g<i/i =1. In general, weak conver- 
gence of measures is not enough to show lower semicontinuity of nonlinear 
integrals or almost everywhere convergence. We claim however that in fact 
the convergence f n — > g takes place strongly in 



lim / \f n (x) - g(x)\dfi(x) = 0. 
Tt - > °° Jm 
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In order to prove this we prove that /„ is a Cauchy sequence in L 1 (d//). We 
take e > and choose N large enough so that 



e 2 



sup\U n (x)-U(x)\ < 



/ U(x)(f n (x) - f m (x))dfJL 
J M 



~ 106 



and 
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S[f n ]<a+- 



hold for n,m> N. Let s(x) = \{f n {x) + f m (x)). Then J M sdfi = 1, s > 0, 
so 

a < £[s]. 

Therefore 



i {£(/„) +S(f m )}-£[s]<^. 



On the other hand, 



j M [/nlog/„ + / ro log/ m ] -slogs} ^ < i TO + 5(/ m )} - + ^ 
so 

1 1 e 2 

2 [/" lo S /n + fm log /m] ~ S log S > d/i < — . 




Denote y = { n . { m and note that — 1 < y < 1 holds « - a.e. Also, elementary 
calculation show that 

^ ^ [/n log fn + fm log / m ] - S log S }> = JC(x) 

holds with 

G( X )=log(l-x 2 ) + xlog^ """ X 



Note that G is even on (-1, 1), that G'( X ) = log = G '(0) = 

and G"{x) = > 2 on (-1, 1). Consequently, 

< X 2 < G( X ) 
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holds for — 1 < x < 1. It follows that we have 

{fn fm) , ^ e 
M Jn ' Jm * 



But, writing I f n — f m \ = y/ f n + fm /} 1? and using the Schwartz inequality 

V Jn~rjm 

we deduce 

/ |/n - fm\dfi < e. 

Therefore the sequence f n is Cauchy in L l (dfi). This proves that the weak 
limit f n d\i — > gd/j, is actually strong /„ — > g in By passing to a 

subsequence if necessary, we may assume also that f n —>g holds also /i- a.e. 
Then, from Fatou's Lemma 



g log gdfx< lim / fjlogfjdp. 
This implies that 

£[g] = a. 

The fact that g solves the Onsager equation follows by taking the Gateaux 
derivative, and thus 

9 = Z-^ bU 

with Z = f M e~ bU dfi. Because U is bounded it follows that g never vanishes 
and because U is Lipschitz, so is g. 

The proof of (B). Let b n — > oo and let us take a subsequence so that g n dfi 
converges weakly to dv. As above, because of ( 1261 1271) and the Arzela-Ascoli 
theorem, we may pass to a subsequence, which we relabel again g n , so that 
U n converge uniformly to a non- negative Lipschitz continuous function U^. 
Let x n be a point where U n {x) attains its minimum U n (x n ) = min xgA f U n (x). 
By passing again to a subsequence we may assume that x n converge to some 
point x. It follows that that U^x) = min m6 jvf U^m) = a. Let be a 
continuous function in M compactly supported in M \ S where S = {m G 
M | Uoolm) = a}. There exists e > so that, for every m in the support of 
(ft, Uooim) > a + 4e. Let us take N so large that sup M \Uoo(m) — U n (m)\ < e 
for n > N and d(x n ,x) < jr Denote a n the minimum of U n . It follows that 
\&n ~ a \ < 2e and U n (m) > a n + e on the support of 0. On the other hand, 
we have 

Z n {b n ) > [ e- bnUn{z) dfx(z) > e~ Kan [ e-^^dfxiz), 
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and using (f2"Tl) we get 

Znipn) > e~ L ce- b "e- bnan 
Therefore, on the support of we have 



(28) 
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and consequently 



n (m) = {Znib))-^-^ 11 ^ < e L c --l e b L ne b n a ne -b n (a n +e) 



(p(rn)g n (rn)djj,(rn) 



M 



<e L c- 1 e b -e- bne \<f>(m)\d(i(m) 

J M 



(29) 



holds, and therefore, as k < 1 we have 



lim / 4>g n d[i = 0. 

n— >oo 1 M 



7 Conclusions 



The generalization of excluded volume interactions, from simple rod-like par- 
ticles to complicated corpora leads to a "simple" equation (the Onsager equa- 
tion) in "complicated" spaces. The examples of single rods, and articulated 
two-rods show significant complexity reduction. The complexity reduction, 
once the problem is phrased correctly, is expected to be generic. The zero 
temperature or high intensity limit of probability distributions of corpora 
concentrates on the minima of certain Lipschitz functions, in general. 



8 Acknowledgments 

Work partially supported by NSF-DMS grant 0504213. I thank A. Zlatos 
and B. Farb for fruitful conversations. 



References 

[1] P. Constantin, Nonlinear Fokker-Planck Navier-Stokes Systems, Com- 
mun. Math. Sci. 3 (2005), 531-544. 



17 



[2] P. Constantin, Smoluchowski Navier Stokes Systems, Contemporary 
Mathematics 429, G-Q Chen, E. Hsu, M. Pinsky editors, AMS, Provi- 
dence (2007), 85 - 109. 

[3] P. Constantin, I. Kevrekidis, E. S. Titi, Remarks on a Smoluchowski 
Equations, Discrete and Continuous Dyn. Syst, 11 (2004) 101-112. 

[4] P. Constantin, I. Kevrekidis, E.S. Titi, Asymptotic states of a Smolu- 
chowski equation, ARMA 174 (2004), 365-384. 

[5] P. Constantin, Ch. Fefferman, E.S. Titi, Global Regularity for Smolu- 
chowski equations coupled to fluids, article in preparation. 

[6] P. Constantin, E.S. Titi, J. Vukadinovic, Dissipativity and GevreyReg- 
ularity of a Smoluchowski Equation, Indiana U. Math J., 54 (2005), 
949-969. 

[7] P. Constantin, J. Vukadinovic, Note on the number of steady states for 
a 2D Smoluchowski equation, Nonlinearity 18 (2005) 441-443. 

[8] P. Constantin, A. Zlatos, The Onsager equation for articulated rods in 
space, in preparation. 

[9] M. Doi, S.F. Edwards, The Theory of Polymer Dynamics, Oxford Uni- 
versity Press, Oxford 1988. 

[10] I. Fatkullin, V. Slastikov, A note on the Onsager model of nematic phase 
transitions, Commun. Math. Sciences. 3 (2005), 21-. 

[11] H. Zhou, H. Wang, M.G. Forest, Q. Wang, A new proof on axisymmetric 
equilibria of a three-dimensional Smoluchowski equation, Nonlinearity 18 
(2005) 2815-2825. 

[12] G. Forest, Q. Wang, R. Zhou, The weak shear phase diagram for nematic 
polymers, Rheologica Acta, 43 (2004), 17-37. 

[13] C. Luo, H. Zhang, P-W. Zhang, The structure of the equilibrium solu- 
tions of one dimensional Doi equation, Nonlinearity 18 (2005), 379-389. 

[14] H. Liu, H. Zhang, P-W Zhang, Axial symmetry and classification of 
stationary solutions of Doi-Onsager Equation on the sphere with Maier- 
Saupe potential. Comm. Math. Sci. 3 (2005), 201-218. 



18 



[15] L. Onsager, The effects of shape on the interaction of colloidal particles, 
Ann. N.Y. Acad. Sci 51 (1949), 627-659. 



[16] F. Otto, A. Tzavaras, Continuity of velocity gradients in suspensions of 
rod-like molecules, SFB Preprint 147 (2004). 

[17] J. Vukadinovic, Inertial Manifolds for a Smoluchowski equation on a 
circle, preprint (2008) 

[18] J. Vukadinovic, Inertial Manifolds for a Smoluchowski equation on the 
unit sphere, preprint (2008). 

[19] A. Zarnescu, Thesis, the University of Chicago, (2006). 



19 



